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ABSTRACT 

We present the first interferometric observations at 18.5 GHz of IRS 1 in NGC 7538. Tliese ob- 
servations inciude images of the nonmetastable ^^NHa (9, 6) masers with a synthesized beam of 2 
arcseconds and images of the continuum emission with a synthesized beam of 150 milliarcseconds. Of 
the maser emission, the previously known feature near i^lsr = — 60kms^^ is spectrally resolved into 
at least two components and we observe several new maser emission features near wlsr — — 57kms~^. 
The new maser emission near — 57kms^^ lies 250 ± 90 milliarcseconds northwest of the maser emis- 
sion near — 60kms~^. All of the masers are angularly unresolved indicating brightness temperatures 
Tg > 2000 K. We are also able to conclusively associate the ammonia masers with the position of 
IRS 1. The excitation of these rare ammonia masers is discussed in the context of the rich maser 
environment of IRS 1 . 

Subject headings: HII regions — ISM: individual (NGC 7538) — ISM: molecules — Masers — Radio 
continuum: ISM — Radio lines: ISM 



1. INTRODUCTION 

The {J,K) = (9,6) -'^^NHs maser was discovered by 
Madden et al. (1986) toward four of 17 Galactic star- 
forming regions surveyed: W51, W49, DR21(0H), and 
NGC 7538. With the exception of confirmation observa- 
tions shortly thereafter, the maser in NGC 7538 has not 
been observed since. The maser in W51 was the sub- 
ject of continued study {e.g., Mauersberger et al. 1987; 
Wilson & Henkel 1988), including the only interferomet- 
ric observations of any these masers prior to the current 
study: a very long baseline experiment by Pratap et al. 
(1991). Empirically, ^''NHa (9,6) masers are known to 
exhibit intensity variability over several months (W49: 
Madden et al. 1986; W51: Wilson, Henkel, & John- 
ston 1990). In the case of the W51 masers, which are 
highly compact, Wilson et al. (1990) and Pratap et al. 
(1991) disagree as to whether the maser gain is satu- 
rated. There is no model for the pump mechanism sup- 
plying the nonthermal population to the (9,6) inversion 
transition for the stimulated amplification of the maser. 
In general, nonmetastable (J > K) transitions are not 
expected to support readily a population inversion, es- 
pecially in the absence of masers in neighboring lines. 
Nevertheless, Brown & Cragg (1991) have studied ortho- 
ammonia {K = 3n, n — 0,1, ...) and found it possible to 
pump an isolated (6, 3) inversion, but they do not extend 
their model to the (9, 6) transition due to computational 
(though not physical) limitations. 

The IRS 1 H II region, assumed by Madden et al. to 
house the ammonia maser in the NGC 7538 complex, 
is particularly rich with maser species, both uncommon 
and ubiquitous (see Pratap et al. 1989; Hoffman et al. 
2003; Galvan-Madrid et al. 2010 for a review). Also, the 
contributions by an accretion disk, jet, and dust to the 
morphology of the continuum emission (see Sandell et 
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al. 2009) and chemistry {e.g., Gibb et al. 2001; Knez et 
al. 2009) are areas of ongoing study (see also Araya et 
al. 2007). In a search for constraints on the complicated 
NGC 7538 region and on the rare masers it hosts, we have 
renewed the study of (9,6) ammonia using the newly 
expanded frequency coverage of the EVLAEI 

2. OBSERVATIONS AND RESULTS 
2.1. 'A' Configuration 

On 2008 October 15 we observed the 18.5-GHz line- free 
continuum emission from IRS 1 using the VLA correlator 
with 600-kHz bandwidth and the 16 antennas outfitted 
with EVLA receivers at the time. The flux and phase 
calibrators used were 3C48 and 2322+509 with approx- 
imately 45 minutes of integration on the IRS 1 target. 
The array was in 'A' configuration resulting in a synthe- 
sized beam of approximately 150 milliarcseconds. The 
primary beam of the VLA antennas at 18.5 GHz is 2.7 
arcminutes. The data were reduced using AIPSH The 
resulting image, shown in Figure [Tl has a background 
noise of 1 m Jy beam" ^ , consistent with instrumental ex- 
pectations. The total imaged fiux density of IRS 1 at 
18.5 GHz is within 10% of 420 mJy. The angular mor- 
phology of the emission is in good agreement with previ- 
ous multi-wavelength studies (see, for example, Sandell 
et al. 2009). 

2.2. 'DnC Array 

On 2010 September 19 we observed the ammonia maser 
line emission (18499.390 MHz rest frequency) using the 

^ The Very Large Array and National Radio Astronomy Obser- 
vatory are facilities of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 

Th e Astronomical Image Processing System is documented at 
|http : //ww w . aips . nrao . edu/ 
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Table 1 

Fitted Spectral Properties of 'SE' 
Maser Feature 



(kms-1) 


I 

(Jybm-l) 


AlIpwHM 
(kms-1) 


-60:21(2) 
-59.79(3) 


0.80(4) 
0.52(6) 


0.49(3) 
0.48(3) 





Note. — The number in paren- 
theses is the uncertainty in the final 
digit. 



EVLA, employing 26 antennas and the WIDAR corre- 
lator. As with the 'A'-configuration observations of the 
continuum emission, the flux, phase, and bandpass cal- 
ibrators used were 3C48 and 2322-1-509 with approxi- 
mately 45 minutes of integration on the IRS 1 target. 
The array was in 'DnC configuration resulting in a syn- 
thesized beam of 4.6 by 1.9 arcseconds at a position an- 
gle of 34 degrees. The observations employed a 500-kHz 
bandwidth centered in frequency on ulsr = — 60kms~^ 
divided into 256 spectral channels resulting in a velocity 
resolution of 0.03 kms~^ and a total velocity coverage 
of approximately 8kms~^. The data were reduced us- 
ing AIPS. The rms background noise in a channel image 
is 17mJybeam~^, consistent with instrumental expecta- 
tions. 

The relative, deconvolved positions of the brightest im- 
age features are precise to within approximately 50 mas. 
The position of the peak maser emission near ?;lsr = 
— 60kms~"'^ lies 250 mas ± 90 mas southeast of the peak 
maser emission near i^lsr = — 57kms^^ and lies 150 mas 
± 70 mas southeast of the location of the peak of contin- 
uum emission from IRS 1. The image of the continuum 
emission from the VLA 'A' observation has been regis- 
tered in absolute position with the image of the contin- 
uum emission from the EVLA 'DnC observation in order 
to summarize these relative positions in Figure [1] The 
absolute position of the images are limited in precision 
to approximately 400 mas by interferometer phase solu- 
tions during data reduction. The velocities in the EVLA 
spectrum corresponding to the different image positions 
are noted in Figure [H 

The masers are not angularly resolved in the EVLA 
observations and have an upper limit on deconvolved 
size of approximately 2 arcseconds. These sizes cor- 
respond to a lower limit on brightness temperature of 
Tb > 2000 K. Assuming that the masers amplify the 
18. 5- GHz background continuum emission from IRS 1 
following Tb = T^gG"'^, we find their gains to be r < —2. 
The spectra of the maser emission features are well sam- 
pled in the EVLA data and show an asymmetric shape to 
the emission near ulsr = — 60kms~^. This shape is well 
fit as the sum of two coincident Gaussians who param- 
eters are listed in Table [T] and whose curves are plotted 
with the data in Figure [2l 

The line emission at the edge of the observing band 
near wlsr — — 57kms~^ has not been detected previ- 
ously (including the 'NW feature in Figure [2]). The ob- 
servations by Madden et al. in 1984 showed no line emis- 
sion at these velocities and so this emission has increased 
in intensity by at least a factor of 50 in the last 27 years. 
We have confirmed this new emission with subsequent 



observations using the Green Bank Telescope that will 
be described by us in a forthcoming paper. This new 
emission also shows asymmetric line shapes indicative of 
several blended narrow maser lines. 

3. DISCUSSION 

In concluding that the ammonia line emission from 
NGC 7538 is maser emission, several assumptions were 
made by Madden et al. (1986) since their arcminute-scale 
angular resolution did not constrain directly either the 
angular extent or the position of the line emission in the 
complex. Namely, it was assumed that the masers were 
associated with IRS 1 and that they were no larger (< 2 
arcseconds) than the molecular clumps observed by other 
authors. With the current data we are able to confirm 
both of these conjectures; definitively associating the line 
emission in position with the IRS 1 continuum emission 
and determining upper limits of approximately 2 arcsec- 
onds for the masering regions. 

In addition, several maser components are spectrally 
resolved, having narrow line widths further suggestive of 
a maser explanation for the emission. The fitted widths 
of these lines are approximately 0.5 km s""'^ which, if cor- 
responding to Doppler shifts from thermal motion in 
the emitting gas, would indicate a kinetic gas temper- 
ature of approximately 90 K. The (9,6) level is approx- 
imately 1100 K above the ground state. The unphysi- 
cally high brightness temperatures of the (9,6) ammonia 
emission lines in IRS 1 compared with their unphysi- 
cally low temperatures calculated from Doppler-width 
assumptions strongly indicate a nonthermal, maser emis- 
sion mechanism for the lines. 

The pumping model of Brown & Cragg (1991), which 
addresses explicitly the (6,3) transition and which is not 
physically inappropriate for the (9,6) transition, requires 
vibrational excitation of the ammonia molecules by in- 
frared blackbody radiation from adjacent dust with a 
temperature of w 300 K. If the kinetic gas temperature 
is also > 300 K then the narrow, sub-Doppler lines widths 
are likely indicative of an unsaturated gain for the maser. 
In contrast, the pumping process suggested by Brown & 
Cragg (1991) readily saturates in terms of the availabil- 
ity of infrared photons. These disparate results leave 
open the following two questions: (1) if the masers are 
unsaturated and if variability is a common hallmark of 
unmet saturation, then why does the intensity of the 
^'LSR — — 60kms~-^ (9,6) maser in NGC 7538 appear to 
have remained unvaried over 30 years? and (2) if, on 
the other hand, the masers are saturated and the gain is 
linear across the line, then why are the line widths so nar- 
row? We note that a similar impasse has been reached in 
the interpretation of the observations of the (9,6) masers 
in W51 by Wilson et al. (1990) and Pratap et al. (1991). 
In the absence of a well tested model for the excitation 
of the (9,6) maser, we compare the current results with 
other maser observations of IRS 1 . 

The angular resolution of the current data permit a 
discussion of the "'^''NHa (9,6) maser in the context of 
the other uncommon H2CO, ^^NHs, and CH3OH masers 
in this region. The recent study of class-II (23.1-GHz) 
CH3OH masers in IRS 1 by Galvan-Madrid et al. (2010) 
shows many similarities with nonmetastable ammonia. 
Both of the main results of the current Letter are also 
seen for the methanol masers: (1) two spectral peaks of 
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Figure 1. The positions of the ^"-NH^ (9, 6) masers in NGC 7538 
IRS 1 from the EVLA 'DnC configuration observations from 2010 
September 19. The crosses mark the positions of two of the maser 
features; the white cross for t^LSR = — 60kms^^ (feature 'SE' in 
Figure[2j and the black cross for i>lsr = — 57kms~^ (feature 'NW 
in Figure [2]l. The sizes of the crosses indicate the uncertainty in 
position of the masers. The locations of CH3OH, ^^NHs, and 
H2CO masers in this region are discussed in §3. The image of the 
continuum emission from IRS 1 is the VLA 'A' configuration image 
at 18.5 GHz from 2008 October 15. The contour levels are -4, 4, 
8, 12, 24, 48, 96, 144, and 176 times the image rms noise level of 
Imjybeam"'^. The beam, plotted in the box at the lower left, is 
171 by 98 milliarcsecond at a position angle of —62 degrees. The 
absolute position uncertainties of the 'DnC' and 'A' observations 
are approximately 400 mas; for display in this figure, the continuum 
peak of the 'A' image has been registered with the continuum peak 
of the 'DnC image from which the maser separations are known. 
The higher precision of the relative positions of the masers from 
each other and from the continuum peak are described in §2. 

maser emission near — 60kms^^ and — 57kms~^ and (2) 
a '^lOO mas northwest-southeast angular separation be- 
tween the peaks. In seeking higher precision than the 
w 100-mas uncertainty in position registration between 
the CH3OH masers and the archival image of the 15-GHz 
continuum emission from IRS 1 used for presentation 
by Galvan-Madrid et al., we have analyzed the archival 
data of the 23.1-GHz observations, including an image 
of the 23.1-GHz continuum emission. Consistent within 
their quoted uncertainties, we find that the northern and 
southern CII3OII masers lie on opposite sides of the peak 
of the 23.1-GIIz continuum emission. The location of the 
more redshifted maser w 40 mas to the north of the con- 
tinuum peak and of the blueshifted maser w 110 mas to 
the south is also seen qualitatively in the current ^''NHa 
distribution, indicating that the two species lie in the 
same kinematic volumes of cloud and outflow, though 
not necessarily in the same clumps of physical conditions. 
Unlike i^'NHa, the 23.1-GHz CH3OH masers show broad 
(2.5kms~^) line widths and little, if any, intensity vari- 
ability over 25 years, likely indicative of saturated maser 



gam. 

The IRS 1 region also exhibits emission from ^^NHs, 
a molecular species that is comparable chemically to the 
^"^NHa presented currently, though likely not similar in 
excitation. Gaume et al. (1991) find a distribution of 
metastable ^^NHs (3,3) 22.8-GIIz masers having many 
similarities with the current presentation of ^"^NHa (9,6) 
masers. In particular, of the numerous ^^NHa masers, 
there are two distinct groupings whose velocities and 
NW-SE position displacement agree well with our Fig- 
ure [1] Two other groupings of ^^NHs masers, for which 
there is no correlated -'^''NHa emission, extend to red- 
der and bluer velocities from which Gaume et al. sug- 
gest a large scale SW-NE velocity pattern of which the 

NII3 velocity distribution presented in our Figure [1] 
is but a small part. We are undertaking wider band- 
width observations of the ^''NHa masers, as well as inter- 
ferometric observations of the ^^NHs masers in a mod- 
ern epoch, that will permit determination of any ad- 
ditional correlations among the new ^''NHa masers at 
'^LSR > — 57kms~^ and the ^^NHs masers which, 22 
years ago, extended over the velocity range —62 kms~^ < 
vlsr < — 52kms~^. Gaume et al. note that the position 
registration of the ^^NHa masers with respect to well- 
modeled images of the continuum emission from IRS 1 is 
plagued by the same absolute position uncertainties we 
have discussed here for other species. Nevertheless, the 
relative location of the ammonia masers and continuum 
emission are well known in the data of Gaume et al. and 
the current data, respectively, and are in agreement with 
each other. 

There are also 4.8-GHz II2CO masers at i'lsr = 
-60.1kms-i and wlsr = -57.9 kms^^ in IRS 1. The 
two interferometric data sets sensitive to both the con- 
tinuum and line emission (Rots et al. 1981; Hoffman et 
al. 2003) find both of the H2CO masers to lie west of 
the 6-cm continuum peak and find the —60.1 kms^^ fea- 
ture to be the westernmost feature. The positions of 
the formaldehyde velocity components stand in contrast 
to the moderate position and velocity correlation among 
CH3OH, i^NHg, and i^NHg in which the -60kms~i fea- 
ture lies to the east of the continuum peak and of the 
redshifted, — 57kms~^ maser. The H2CO masers are 
variable: there is long-term intensity variability in the 
H2CO masers which may be correlated with the appear- 
ance of the new ^^NHa masers reported in this Letter. 

Common temporal variability of various maser species 
may provide insight into their various excitation condi- 
tions. For example, two maser species may exist at dif- 
ferent densities or velocities or positions, though both 
may be dependent on the same infrared source for radia- 
tive excitation and therefore exhibit common variabil- 
ity. Or, the relatively slow passage of a bulk material 
disturbance through the maser environments will affect 
only those masers sufficiently sensitive to changes in den- 
sity or collisions when the disturbance arrives at their 
location. We examine possible excitation scenarios in 
IRS 1 in order to explain two differing aspects of the 
^^NHa masers reported here: (1) the unchanged inten- 
sity of the component near ulsr = — 60kms~"'^ and (2) 
the > 50-fold increase in intensity of the component near 
"LSR = -57kms"^ 

Neither the CH3OH nor ^^NHs masers discussed above 
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Figure 2. Spectrum of the ^^NHs (9, 6) masers in NGC 7538 IRS 1 from the EVLA 'DnC configuration observations from 2010 September 
19 (solid histogram) after subtraction in tiie u,v plane of the contribution from continuum emission. The 'SE' and 'NW labels correspond 
to the deconvolved locations of the emission in the image plane of Figure [T] This single spectrum includes all of the maser emission since 
all of the emission is convolved within the synthesized beam on the sky of size 4.6 by 1.9 arcseconds at a position angle of 34 degrees. The 
dotted and dashed lines are the two best-fit Gaussians whose sum is the asymmetric line profile near tiLSR = — 60kms~^, as summarized 
in Table [Tl and described in §2. 



have exhibited any intensity variabiHty in either the 
— 60kms~^ or — 57kms~^ velocity components since 
their discovery approximately 30 years ago (Wilson et 
al. 1984; Johnston et al. 1989). Since class-II methanol 
masers are thought to depend on excitation from infrared 
radiation, the current ^"^NHa variability results are not 
suggestive of a common radiative excitation for ammo- 
nia and methanol. Hoffman et al. (2003) note common, 
long-term variability of H2CO and H2O masers, further 
described by Lekht et al. (2003, 2004; see also Hoffman 
et al. 2007). Araya et al. (2007) have continued to moni- 
tor the variability of the two H2CO velocity components, 
suggesting possible successive excitation by a moving dis- 
turbance or jet similar to that described by Araya et 
al. (2009) to explain the behavior of CH3OH masers in 
DR21(0H). It is possible that the ^^^NHs maser near 
'^LSR — — 57kms~-'^ has been excited sometime in the 
past 27 years by the arrival of a bulk material distur- 
bance. Such a disturbance, if responsible for all of the 
long-term variability in H2CO, NH3, and H2O, would 



need to be an isotropic flow or jet directed both north 
and south from the center of IRS 1 since the H2CO maser 
and the variable nonmetastable ammonia maser lie on 
opposite sides of the peak of continuum emission. Ongo- 
ing monitoring observations are essential for determining 
the presence of an excitation disturbance traversing the 
region. Since Brown & Cragg (1991) find the pumping 
of the (6,3) line, and likely also the (9,6) line, to require 
both radiative and coUisional excitation, the existence of 
the rare ammonia, formaldehyde, and methanol masers 
in IRS 1 may be indicative of a relatively uncommon 
superposition of infrared radiation and material outflow. 

The radio continuum morphology of IRS 1 is extended 
north-south and has been successfully modeled as a bipo- 
lar outflow with the southern half more blueshifted than 
the north (for example, Sandell et al. 2009). The vari- 
ous maser species discussed above have been speculated 
inconclusively to include some or all of a jet, a shock, 
microwave free-free emission, infrared dust emission, and 
molecular clumps for their excitations. Since relative po- 
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sitions within a data set are generally more precise than 
comparisons of absolute positions between data sets, we 
have restricted our discussion in this Letter to the rel- 
ative positions of the NH3, H2CO, and CH3OH masers 
and C- and K-band images of the continuum emission. 
Nevertheless, in addition to the broad correlations we 
have noted, other authors have suggested several spe- 
cific associations among individual masers (additionally 
H2O and OH) and among thermal and continuum fea- 
tures that are supported by precise absolute positions 
{e.g., Hoffman et al. 2003; Galvan-Madrid et al. 2010). 
Wc plan to explore possible associations and interrelated 
physical conditions in a forthcoming paper. 

4. CONCLUSION 

We report the discovery of new maser emission from 

the nonmetastable (J, K) = (9, 6) inversion transition of 
"NH3 in NGC 7538 IRS 1. Both the previously known 
feature near t^LSR = ^60 kms^^ and the new masers near 
— 57kms~^ are imaged interferometrically for the first 
time, showing different positions in IRS 1 in correlation 
with some other maser species. The brightness tempera- 
tures {Tb > 2000 K) and line widths (O-Skms^i) of the 
features demonstrate conclusively the nonthermal nature 
of the emission. We plan future EVLA observations of 
these masers with improved angular resolution in order 
to constrain their excitation mechanism. 

This work is supported by the Thomas Penrose Ben- 
nett Prize Fund and the Lovejoy Science Fund of St. 
Paul's School. We arc grateful to the NRAO staff in 
Socorro for flexible scheduling and guidance during the 
EVLA transition. We are indebted to an anonymous 
referee for comments that significantly improved the 
manuscript. 
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